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The discovery of transiting circumbinary planets by the Kepler mission suggests that plan¬ 
ets can form efficiently around binary stars. None of the stellar binaries currently known 
to host planets has a period shorter than 7 days, despite the large number of eclipsing bi¬ 
naries found in the Kepler target list with periods shorter than a few days. These compact 
binaries are believed to have evolved from wider orbits into their current configurations 
via the so-called Lidov-Kozai migration mechanism, in which gravitational perturbations 
from a distant tertiary companion induce large-amplitude eccentricity oscillations in the 
binary, followed by orbital decay and circularization due to tidal dissipation in the stars. 
Here we explore the orbital evolution of planets around binaries undergoing orbital decay 
by this mechanism. We show that planets may survive and become misaligned from their 
host binary, or may develop erratic behavior in eccentricity, resulting in their consumption 
by the stars or ejection from the system as the binary decays. Our results suggest that 
circumbinary planets around compact binaries could still exist, and we offer predictions 
as to what their orbital configurations should be like. 

planet dynamics and stability — close stellar binaries — multiple stellar systems 


To date, the Kepler spacecraft has discovered eight binary star sys¬ 


terns harboring ten transiting circumbinary planets Doyle et al. 

201 Ij; 

Welsh et al.|(|2012j;|Orosz et al. 

(|2012|?J; Schwamb et al. 

2013j; 

Kostov et al.|(|20131|2014J;|Wels 

1 et al.|(|2014|(. These sys- 


tems have binary periods ranging from 7.5 to ~ 41 days, while the 
planet periods range from ~ 50 to ~ 250 days. Remarkably, no 
transiting planets have been found around more compact stellar bi¬ 
naries, those with orbital periods of < 5 days. Planets around such 
compact binaries, if orbiting in near coplanarity, should have tran¬ 
sited several times over the lifetime of the Kepler mission. How¬ 
ever, the shortest period binary hosting a planet is Kepler-47(AB) 
with 7.44 days, despite the fact that nearly 50% of the eclipsing 
binaries in the early quarters of Kepler data have periods shorter 
than 3 days |Slawson et al.| j2011j . Thus, the apparent absence 
of planets around short-period binaries is statistically significant 
[e.g. [Armstrong et^ ( |2014j ]. 

It is widely believed that short-period binaries (< 5 days) 
are not primordial, but have evolved from a wider configurations 


via Lidov-Kozai (LK) cyclesILidov 

l(fT962 1 ; Kozai|(11962 1 with tidal 

friction 

Mazeh & Shaham 11979J; 

Eggleton & Kiseleva-Eggleton| 

pO^; 

Fabrycky & Tremaine J2007l. This “LK-Hide” mechanism 


requires an external tertiary companion at high inclination to ex¬ 
cite the inner binary eccentricity such that tidal dissipation becomes 
important at pericenter, eventually leading to orbital decay and cir¬ 
cularization. A rough transition at an orbital period of 6 days has 
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been identified as the separation between “primordial” and “tidally 
evolved” binaries [Fabrycky & Tremainej ( |2007| >. Indeed, binaries 
with periods shorter than this threshold are known to have very 
high tertiary companion fractions [of up to 96% for periods < 3 
days; see [Tokovinin et dl]j2006j ], supporting the idea that three- 
body interactions have played a major role in their formation. 

In synthetic population studies jFabrycky & Tremainej 
l |2007j , stellar binaries with periods shorter than ~ 5 days evolved 
from binaries with original periods of ~ 100 days. Interestingly, it 
is around binaries with periods < 100 days that transiting planets 
have been detected. It is thus plausible that current compact bina¬ 
ries with a tertiary companion may have once been primordial hosts 
to planets like those detected by Kepler. 

In this work, we study the evolution and survival of plan¬ 
ets around stellar binaries undergoing orbital shrinkage via the 
“LK-l-tide” mechanism. We follow the secular evolution of the 
planet until binary circularization is reached and binary separation 
is shrunk by an order of magnitude. We show that the tertiary com¬ 
panion can play a major role in misaligning and/or destabilizing the 
planet as the binary shrinks. 

A PLANET INSIDE A STELLAR TRIPLE 

Consider a planet orbiting a circular stellar binary of total mass 
Min = rno 4- mi and semimajor axis ain; the binary is a member 
of a hierarchical triple, in which the binary and an outer companion 
of mass Mont orbit each other with a semimajor axis ttout fflin. 
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Steady-state planetary orientation for a eireular, preeessing inner binary 
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b) 


c) 
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Figure 1. The planet’s angular momentum orientation Lp processes around the equilibrium solution Lp^eq (red vector), which is obtained from balancing 
the torques acting on the planet due to the inner binary and the outer companion. The limiting cases are (a) Lp^eq || Lm (planet at small distance Up from 
the inner binary, where outer torque is negligible) and (c) Lp,eq || Lout (at large distance, where the inner torque is negligible), while (b) represents the 
intermediate cases where the two torques have comparable magnitudes. The continuous variation of Lp^eq with Up (always coplanar with Lm and Lout) 
going from (a) to (b) to (c) defines the so-called “Laplace surface”. When the inner binary axis Lin precesses slowly around Lout by the action of the outer 
companion, Lp^gq follows adiabatically and stays coplanar with Lin and Lout- If the torque from the inner binary is slowly decreased in time, (e.g., due to 
orbital decay), Lp,eq not only will process around Lout, but will also change its inclination, going smoothly from regime (a) to regime (c). 


The secular (long term) gravitational perturbations exerted on the 
planetary orbit from the quadrupole potential associate(0with the 
inner binary and that from the outer companion cause the two vec¬ 
tors that determine the orbital properties of the planet, the angular 
momentum direction Lp and the eccentricity vector Op, to evolve 
in time. The inner binary tends to make Lp precess around Lin, 
the unit vector along the inner binary’s angular momentum, at a rate 
approximately given by 




( 1 ) 


where Op is the semi-major axis of the planet, rip = 
is the planet’s mean motion frequency (assumed to be on a circular 
orbit), and /iin = momi/Min is the reduced mass of the inner 
stellar pair. Similarly, the outer companion of mass Mout tends to 
make Lp precess around Lout at a rate approximately given byq 




p-out 


= rip 


Moot 

Min 


/ flp 
\ dout 
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( 2 ) 


In general, when the torques from the inner binary and the outer 
companion are of comparable magnitude, Lp will precess around 
an intermediate vector Lp,eq, which corresponds to the equilib¬ 
rium solution (i.e., d Ijpjdt = 0) of the planet’s orbit under the 
two torques. For a general mutual inclination angle iin-out between 
the inner and outer orbits (where cosiin-out = Lin • Lout), the 
equilibrium inclination of the planet (the so-called “Laplace sur¬ 
face”; see [Tremaine et al,] ( |2009[ >; [Tamayo et al.| ( |2013) ), can be 


^ If the inner binary has an equal-mass ratio and the outer companion has 
zero eccentricity, the octupole-order terms in the potential vanish exactly. 

^ Although we assume a circular outer companion here, the eccentricity 
of the outer orbit Bout can be taken into account by replace Uout with 

Uout^l-e^^t . 


found as a function of its semimajor axis, for which Lp,eq is al¬ 
ways coplanar with Lin and Lout, with limiting states correspond¬ 
ing to alignment with the inner binary (i.e., Lp,eq || Lin) at small 
Op, and alignment with the outer companion (i.e., Lp_eq |j Lout) 
at large Op. The transition between these two orientations happens 
rapidly at the so-called “Laplace radius” vl, obtained by setting 
fip.out = f2p,in, and is given by 


tl 


f^in 

2A(fout 



1/5 


(3) 


Fig. [T] illustrates the three regimes of the planet’s equilibrium ori¬ 
entation: (a) flp-in flp-out (binary-dominated regime, or Op <C 
tl); (b) flp-out ~ flp-in (transition regime, or ttp ~ rr,); and (c) 
flp-out S> flp-in (companion-dominated regime, or Up ^ vl)- 

In general, however, the vector Lin is not fixed in space, but 
slowly precesses around Lout, owing to the torque from the outer 
companioij^ This means that the plane normal to Lin x Lout, 
where the equilibrium orientation vector Lp,eq lives, is slowly ro¬ 
tating (Fig.[^. This rotation rate is of order 


Din-out — rt-in 


Mo- 


Mi, 



(4) 


where riin = a\^ is the mean motion of the in¬ 

ner binary. Note that Dp-out/Din-out = (flp/ain)^^^ » 
1 in the companion-dominated regime, and Dp-in/Din-out = 
(Dp-in/Dp-out)(ap/ain)^^^ = (rL/ap)'^(ap/oin)®^^ » 1 m the 
binary-dominated regime. This means that the precession of Lin is 
always slow enough for Lp to adiabatically follow. In other words, 
the classical Laplace equilibrium formalism remains valid in the 


^ Strictly speaking, both Lin and Lout precess around the total angular 
momentum vector of the system; however, for the hierarchical configura¬ 
tions presented here, the outer orbit contains most of the angular momentum 
of the system, implying that Lout is approximately fixed in space. 
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Figure 2. The three relevant precession frequencies (flp.in, Op-out and Oin-out) as a function of the shrinking binai’y semimajor axis ain. The binary starts at 
semimajor axis ajn q = 0.3 AU and circularizes at Qin,/ = 0.024 AU (vertical black line). The other parameters are = Mout = IMq, /im = 0.25, 
(lout = 30 AU and Cout = 0. The left panel shows the case with Op = 2 AU, and the right panel one with Op = 1 AU. On the left panel, ain crosses 
“in,L = 0.097 AU (thick vertical gray line), during orbital decay, then the planet transitions from the binary-dominated regime into the companion-dominated 
regime. On the right panel, ctin,/ > ciin,L = 0.017 AU) and the planet will stay in the binai'y-dominated regime throughout the binary orbital decay. Note, 
that in this example, Up = 1 AU is very close to the initial binary, and dynamical instabilities (not captured by secular calculations) might make the survival 
of these planets difficult during the early Lidov-Kozai cycles of the binary. 


frame corotating with Lin, and the three vectors Lin, Lout and 
Lp,eq remain coplanar at all times. Since the evolution is adiabatic, 
if Lp starts parallel to Lp,eq, it will remain parallel to the evolving 
Lp^eq at later times, provided that this equilibrium orientation is a 
stable solution |Tremaine et al.| ( [200^ . 

As studied by Tremaine et al.| ( [20Q9) , when fin-out > 69°, 
circular orbits on the Laplace surface are unstable to linear per¬ 
turbations in the planet’s eccentricity vector Op vector for a range 
of Up around vl- This instability manifests itself as an exponen¬ 
tial growth of Cp, until non-linear effects come into play, resulting 
in erratic behavior in both inclination and eccentricity. This means 
that above this critical value of fin-out, planets cannot be placed at 
Up ~ fi, since the resulting high eccentricities could bring them 
too close to the binary, at which point they may collide with the 
central stars or be ejected from the system [e.g.,|Holman & Wiegert| 
( [19^ ]. 


Now consider what will happen to the planet’s orbit as the 
inner binary undergoes orbital decay. For simplicity, let us assume 
that the binary remains circular during this process, and that the an¬ 
gle fin-out remains unchanged. Since orbital decay takes place over 
a time scale fdecay much longer than the other relevant time scales 
(l/^p -ill) 1/^p-out and l/riin-out), the system will evolve adiabat- 
ically. Thus, if the planet initially resides in the binary-dominated 
regime (flp-in ^ flp-out, a-p <C ^l), and lives on the Laplace sur¬ 
face ( Lp II Lp^eq), it will transition to the companion-dominated 
regime (flp-in ^ flp-out, a-p ^ fi) through the intermediate stage 
(flp-in ~ flp-out), as the inner binary’s semi-major axis de¬ 
creases. For a given value of Up, the transition occurs when ai„ 
passes through a critical (“Laplace”) value. 


ain.L = 0.017 


fMo 


\ 4)/ij 


1/2 


1 AU 
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^^■out 

30 AU 


- 3/2 


AU 


(5) 

obtained by replacing ri, = Up in Eq.|^and solving for ain. If the 
transition region (Op ~ tl) is stable, we expect the planet’s orbit to 
evolve smoothly following the Laplace surface [i.e., (a)—^-jb)—^(c) 
in Fig.[T). For fin-out > 69°, however, the planet will encounter 


an instability when Up ~ tl, and may undergo erratic evolution, 
which may result in the planet being destroyed or ejected. 

In the “LK-l-tide” scenario for the formation of compact bi¬ 
naries, the final inner binary separation flin,/ depends on the prop¬ 
erties of the outer companion (Mout, flout and the initial inclina¬ 
tion fin-out) as well as on the short-range force effects between the 
inner binary members |Fabrycky & Tremaine| pOOV^ . Thus, for a 
given stellar triple configuration, the inner binary may or may not 
reach down to flin,L, depending on the value of Op (see Fig. |^. If 
flin,/ > flin.L, Or equivaleutly, if 


flp < 1.26 


/Mo„t\ 

\ 4/tin / 


\ 3/5 
Clout \ 

30 AU J 


0-in,f \ 

0.03 AU J 


2/5 


AU , 


( 6 ) 

the planet will never cross the intermediate regime (Op ~ tl), 
and it will thus remain “safe” (stable), regardless of the inclination 
fin-out , surviving the orbital decay of its host binary. 


EVOLUTION OF PLANETARY ORBITS AROUND 
BINARIES UNDERGOING LIDOV-KOZAI CYCLES WITH 
TIDAL FRICTION 


The greatest caveat to the application of classical Laplace equilib¬ 
rium is that the inner binary does not remain circular during orbital 
decay. Indeed, in the “LK-l-tide” mechanism Eggleton & Kiseleva-| 
|Eggleton] ( |2001| l; |Fabrycky & TremaineH2007^ the inner binary ex- 
hibits large oscillations in inclination and eccentricity under the in¬ 
fluence of the external stellar companion. Thus the binary axis Lin 
not only processes around Lout, but also undergoes nutation. The 
variation of the inner binary’s eccentricity vector ein also affects 
the torque on the circumbinary planet. 

To track the evolution of the planet’s orbit during the LK 
oscillations and orbital decay of the inner binary, we solve numer¬ 
ically the secular equations of|^f the planet’s eccentricity vector 


The secular equations of motion govern the evolution of the orbital ele¬ 
ments instead of the position and velocity of individual bodies. 
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Figure 3. An example of the coupled evolution of an inner binary within a 
stellai* triple (mo = mi = O.SMq, Mout = IMq, aout = 18 AU and 
initial ain,o = 0-3 AU and iin-out,0 = 73°) plus a planet with semimajor 
axis ttp = 1.5 AU. The different panels show: (top) binary semimajor axis 
ain (green) and planet semimajor axis ap (blue) and the Laplace radius 
(cyan); (top middle) eccentricity of the binaiy ein (green) and eccentric¬ 
ity of the planer ep (blue); (top bottom) inclination of the binary iin-out 
(green) and inclination of the planet Zp-out (blue) with respect to the outer 
companion; and (bottom) mutual inclination between the planet and the bi¬ 
nary Zp-in- The eccentricity and inclination of the binary exhibit LK cycles 
for about 10% of the integration time, until short-range forces arrest these 
oscillations (freezing ein at high values), at which point a slow phase of 
orbital decay takes place. The planet starts in the binary-dominated regime 
(czp/^L,o = 0.65). Its inclination Zp-out follows closely that of the inner 
binary zin-out until crosses ap (note that the definition of in Eq. 
does not take into account the eccentricity of the inner binary ein), at which 
point these two inclination angles decouple from each other. The planet ends 
in the companion dominated regime {aplr^j = 1.3), and its inclination 
with respect to the binary Zp-in eventually settles into a constant value of 
32°. 

Gp and angular momentum vector axis Lp (see Supplementary In¬ 
formation), along with the evolution equations of the stellar triple. 
We use the formalism of |Eggleton et al.|jl998| t to follow the in¬ 
ner binary’s orbit and parametrize the stellar tidal dissipation rate 
using the weak friction model with constant tidal lag time. In the 
following, we focus on a few representative examples and discuss 
the general behavior for the evolution of the four-body system. 

Fig. [^depicts a system where the stellar triple has param¬ 
eters mo = mi = 0.5M©, Mout = IM©, Uout = 18 AU and 
eout = 0 and initial values ain,o = 0.3 AU and Zin-out,o = 73°, and 
where the planet is initialized on a circular orbit at Up = 1.5 AU 
with Lp aligned with Lin. The parameters for the inner binary 
and the planet are chosen to roughly correspond to the discov¬ 
ered Kepler systems. The parameters for the outer orbit are cho¬ 
sen to ensure that LK cycles are not suppressed by short-range 
forces and to guarantee the efficient orbital decay of the inner 



[AU] 


Figure 4. Classical Laplace equilibrium surface (valid for ein = 0) at 
the beginning (red curve) and after circularization of the inner binary (or¬ 
ange curve) for the triple configuration of Fig. (mo = = O.SMq, 

Mout = 1M0, ain_o = 0.3 AU, iin-out,o = 73°, = 0.053 AU, 

iin-out,/ = 46.1° and Oout = 18 AU). The dotted portion of the red 
line indicates the range of radii at which the equilibrium surface is unsta¬ 
ble ^Tremaine et al. [20091 T’he final Laplace surface is stable for all ap 
since ip-out,/ < 60°. The vertical dash-dotted lines indicate the Laplace 
radii at the beginning (rj^ o) and end of the binary orbital evolu¬ 

tion. For the different values of ap, vertical an'ows connect the initial and 
final states, representing the evolution of the planet’s inclination obtained 
from the numerical calculations. In each case, the planet orientation is ini¬ 
tially aligned with the local Laplace surface (or approximately aligned with 
the binary for ap < 1.5 AU); after the inner binary has decayed and cir¬ 
cularized, the planet inclination settles into a value coincident with final 
Laplace surface. Note that, for illustrative purposes, we include values of 
ap down to 1 AU; however, dynamical stability dictates that only planets 
outside ap ~ dam ~ 1.2 AU (when ein ~ l; |Flolman & Wiegert|1999) 
should survive the early LK cycles of the inner binary. 

binary [Fabry cky & Tremain^ ( |2007) . In our calculations, the oc- 
tupole term in the potential has been ignored in the evolution equa¬ 
tions of the planet and the inner binary, a justified simplification 
since mo = mi and Bout = 0. The inner binary experiences LK 
oscillations and circularizes within a Hubble time provided that 
enough tidal dissipation is present in the stars. The final (circu¬ 
larization) semimajor axis is atn,/ = 0.053 AU (corresponding 
to an orbital period of 4.5 days). In this example, the planet ini¬ 
tially resides in the binary-dominated regime, with flp-in/Hp-out ~ 
65 (flp/AU)”® « 8.6, and ap/rL,o ~ 0.65. After the inner bi¬ 
nary has circularized, the planet lies in the companion-dominated 
regime, with flp-in/flp-out ~ 0.27 and ap/r_L j « 1.3. We see that 
the planet remains on a circular orbit throughout its entire evolu¬ 
tion, despite the large variations in ein during the LK cycles. The 
longitude of nodes of the planet (not shown in the figure) closely 
follows that of the inner binary during the early LK cycles and after 
circularization, implying that for a large fraction of the time Lp is 
coplanar with Lin and Lout- The planet’s inclination ip-out also 
follows the binary inclination iin-out during the early stage of the 
LK cycles (third panel from top), but it decouples from the inner 
binary after ain has started decreasing. At the end of the integra¬ 
tion, when the binary has circularized, the binary and planet are 
misaligned by 32° (fourth panel) and the planet inclination has set¬ 
tled onto a steady-state value. This final value ip-out — 14° agrees 
with the equilibrium value of the end-state Laplace surface (with 
fflin,/ = 0.053 AU and im-out,/ = 46°) evaluated at Op = 1.5 AU. 

We have carried out calculations for a range of values of 
Up for the same stellar triple configuration of Fig. [^ The results 



















































Survival of Planets Around Shrinking Stellar Binaries 5 



10 “ 10 ^ 10 ® 10 ^ 10 “ 10 ® 10 “ 10 ^ 10 ^ 10 ' 10 “ 10 ^ 

Time [yr] Time [yr] 


Figure 5. Similar to Fig.j^ but for a triple system that is initialized at a higher inclination *in-out,0 = 83°. The other parameters for the tare ain o = 0-3 AU 
and Uout = 30 AU, with the same stellar masses as in Fig. Two examples are shown: ap = 1.2 AU (left panels) and ap = 1.8 AU (right panels), 
both exhibiting quite different planetary evolution compared to Fig.|^ In the case of ap = 1.2 AU, the planet starts with apjri^fi = 0.38 and ends with 
apIvLj = 1.02; in the case of ap = 1.8 AU, the planet starts with apjri^Q = 0.57 and ends with apjrj^f = 1.54. 


of these calculations are summarized in Fig. which shows the 
Laplace surfaces at the beginning and at the end of the evolution, 
when the inner binary is circular. In each case, the planet is ini¬ 
tially aligned with the equilibrium orientation Lp,eq, which is in 
near alignment with the inner binary for ap < 1.5 AU. We find 
that the planet’s inclination evolves smoothly for all these cases as 
the binary experiences LK oscillations and orbital decay. Despite 
the complexity of the “intermediate” states, in which the binary 
develops large eccentricities and the standard Laplace equilibrium 
is not well defined, we find that in the end, the planet’s inclina¬ 
tion always lands on the final Laplace surface. Thus, these planets 
survive the orbital decay of the inner binary, but become inclined 
respect to it by an angle given by ip-in,/ = hn-out,/ - ip-out, f (with 
*in-out./ ~ 46° for the parameters adopted in Figs.[^and[^, where 
ip-out,/ matches the equilibrium inclination of the final Laplace 
surface. Since the Laplace equilibrium inclination angle decreases 
with increasing ap, we predict that the angle ip-m of the planets that 
survive will increase monotonically with increasing ap. 

As noted before, when the mutual inclination iin-out be¬ 
tween the inner circular binary and the external companion is 
greater than 69°, a portion of the Laplace surface is unstable 
[Tremaine et al.| ( |2009^ . In principle, a circumbinary planet may 
suffer a similar instability as a binary with large initial fin-out un¬ 
dergoes “LK-l-tide” orbital decay. In Fig.[^ we show two examples 
(ap = 1.2 AU and ap = 1.8 AU for the left and right panels, re¬ 
spectively) of planets within a stellar triple with flout = 30 AU, 
flin.o = 0.3 AU, and fin-out,o = 83° (the other parameters are the 
same as in Fig.[^. At this initial inclination, the inner binary attains 
very high eccentricities and can circularize very efficiently (alter¬ 


natively, it requires relatively small tidal dissipation in the stars to 
circularize within a Hubble time). The final binary separation is 
flin,/ = 2.55 X 10“^ AU (period of 1.5 d) and the inclination an¬ 
gle freezes out at fin-out,/ = 65.3°. The behavior of the planets 
is markedly different from the one depicted in Fig.[^ For a planet 
located at flp = 1.2 AU (left panels of Fig.|^, the inclination an¬ 
gle fp-out does not evolve smoothly as the inner binary decays, but 
suffers a jump as vl crosses flp, subsequently oscillating around a 
reference angle. Moreover, the orbital eccentricity rapidly grows 
until it starts oscillating around a mean value of (cp) ~ 0.16, 
maintaining from then on a steady-state behavior. For a planet at 
flp = 1.8 AU (right panels), the orbital evolution is even more 
complex. In this case, the exponential growth in eccentricity does 
not saturate at a moderate value. Instead, Cp reaches values close to 
1. The erratic evolution in Cp is accompanied by a similar behav¬ 
ior in the planet’s inclination fp-out - Instead of oscillating around a 
mean (equilibrium) value, fp-out covers the entire range (0°, 180°). 
The high planet eccentricities reached in this case make it very un¬ 
likely for the planet to survive the orbital decay of the inner binary. 
Such high eccentricities will inevitably bring the planet too close 
to the inner binary, a region that is known to be unstable |Holman| 
|& Wieg^ jl999| l; [Mudryk & Wu| ( |2006^ . In this case, ejections 
from the system or physical collisions with the central stars are to 
be expected. 

In Fig. 1^ we show the initial and final inclinations (top 
panel), and the respective final eccentricities (bottom panel), com¬ 
puted for a set of values of flp using the same stellar triple con¬ 
figuration of Fig. 1^ As in Fig. these results are shown to¬ 
gether with the Laplace equilibrium surface solutions for the ini- 
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Figure 6. Top panel: similai' to Fig. but for the triple configuration 
shown in Fig.|^(ain^o = 0.3 AU, iin-out,0 = 83°, tiin./ = 0.026 AU, 
*in-out,/ = 65.3°). Bottom panel: the correspondingeccentricity of 
each planet. Error bars specify the oscillation amplitude around the mean 
value. Purple symbols denote planet orbits that were excited into eccentric 
states. Wide blue bands (for Op > 1.5 AU) denote planet orbits with erratic 
behavior in inclination, covering the entire range [0°, 180°), at any point 
during their evolution. Note that, for some values of Up, planets do reach 
regular values in eccentricity and inclination even after having experienced 
erratic evolution during a finite period of time before binary circularization; 
such cases are still depicted by blue bands, since their survival is deemed 
unlikely (see Fig S2). 

tial state (ain.o = 0.3 AU and iin-out,o = 83°) and the final state 
(din,/ = 0.026 AU and iin-out,/ = 65.3°). Unlike Fig.[^ we find 
that, depending on Op, planet orhits do not always stay circular, 
and their inclinations ip-out do not always land exactly on the fi¬ 
nal Laplace surface. For Op < tlj, planets end up very close (on 
average) to the final Laplace surface (while exhibiting some mi¬ 
nor oscillations around it), and maintain a negligible eccentricity. 
For Up > tlj, planets suffer a small kick in eccentricity as they 
cross the “transition” regime (Up = vl), and their inclinations os¬ 
cillate with significant amplitude around a mean value that is close, 
but not necessarily equal to, the one given by Laplace equilibrium 
(see Fig. left panel). At even larger Up (> 1.5 AU), we find 
that the evolution of the planet is no longer regular (see Fig. 
right panel): both Cp and ip-out undergo large-amplitude, erratic 
variations (gp ~ 0 - 1 and Ip-out — 0°-180°). Indeed, for large 
values of Up, the planet’s evolution is most likely chaotic, since 
the results depend sensitively on the initial conditions (see Supple¬ 
mentary Information and Fig. S2). Erratic evolution (in eccentricity 
and inclination) may last indefinitely or may end before circular¬ 
ization of the inner binary has completed, in which case planets 
can exit the erratic phase at a random inclination (including angles 
> 90°). In either case, these planets, having experienced erratic, 
large-amplitude variations of Cp, are likely to be ejected from the 
system or to collide with the binary stars. 

In the above, our calculations have ignored the mass of the 
circumbinary planet nip based on the assumption that the plan¬ 


etary mass is always much smaller than Min and Mout- How¬ 
ever, over secular time scales, a finite planet mass can affect the 
dynamics of the inner binary to the point of suppressing the ec¬ 
centricity oscillations caused by the tertiary [Holman et al.| ( |1997) l. 
The planet-induced precession frequency of the binary is of order 
flin-p — riin {vTlpj Aljn) (Uin/Up) . The Condition Din-p — Hin-out 
allows US to define the critical planet mass 


mp,crit 0.13Mj 


(MouA 

( ap Y 

/ Qiout \ 

[iMqJ 

V1.5AU/ 

V30AU/ 


(7) 


(where Mj is the mass of Jupiter) above which the precession of 
the binary due to the planet is faster that due to the tertiary star. 
For small mp, the effects of a finite planet mass on the LK cy¬ 
cles are qualitatively similar to those of other short range forces 
[Fabrycky & Tremaine| ( |2007t ; |Liu et al.| ( |20r5t , imposing an up¬ 
per limit on the maximum eccentricity of the binary. In Fig. [^ we 
show the maximum eccentricity achieved by the inner binary as a 
function of planet mass obtained from integrations of the 4-body 
secular system (see Supporting materials). For the example de¬ 
picted in Fig. [^(Up = 1.5 AU and mp,crit — 0.6Mj), we find 
that nip > IMj ~ 1.7mp,crit is enough to substantially suppress 
the oscillations in ein. For nip < 0.3Mj ~ 0.5mp,crit (about the 
mass of Saturn), the minimum pericenter separation of the binary 
ain,o(l 6in,inax ) « 0.09ain,o is only 17% larger than 0.077ain,o, 
the value corresponding to nip = 0. Such a planet {nip < 0.3Mj) 
will only delay the orbital shrinkage of the inner binary, but not 
prevent it (see Supplementary materials for an example). 

Throughout this paper, we have included only the 
quadrupole potential from the tertiary companion acting on the in¬ 
ner binary and the planet. This is a good approximation when the 
companion has zero orbital eccentricity. For general companion 
eccentricities, octupole and higher-order potentials may introduce 
more complex dynamical behaviors for the inner binary and for the 
planet (see, e.g.,|Ford et al.|(|2000|(;|Naoz et al.|(|2011|(;|Katz et al.j 
p0TT) ;[L 


.lu et a 


1 2015 l). For example, in H-body calculations 
(which include high order terms automatically) the planet may at¬ 
tain a non-zero eccentricity as the inner binary decays even in the 
moderate inclination case (see the Supplementary material for one 
such example). A systematic study of these complex “high-order” 
effects is beyond the scope of this paper and will be the subject of 
future work. 


DISCUSSION 

We have explored the orbital evolution of planets around binaries 
undergoing orbital decay via the “LK-l-Tide” mechanism driven by 
distant tertiary companions. We have shown that planets may sur¬ 
vive the orbital decay of the binary for tertiary companions at mod¬ 
erate initial inclinations (iin-out,o ^ 75°). In such case, planets 
on circular orbits adiabatically follow an equilibrium solution as 
the triple system evolves, becoming misaligned with their host bi¬ 
nary; the final misalignment angle ip.in is a monotonically increas¬ 
ing function of the binary-planet distance ap. At higher inclinations 
(*in-out,o ^ 80°), the adiabatic evolution is broken when plan¬ 
ets encounter an unstable equilibrium. Then the planet orbit can 
develop erratic behavior in eccentricity and inclination. Very ec¬ 
centric circumbinary orbits may be disrupted by the inner binary 
via dynamical instabilities, resulting in either the ejection of the 
planet or its collision onto the stars. Interestingly, even in this high- 
inclination regime, we have found that some planets may evolve 
into stable, misaligned and eccentric orbits. 
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Figure 7. Left: maximum eccentricity ein.max of the inner binary in the triple configuration of Fig.j^achieved during the LK cycles as a function of planet 
mass mp for three different values of the planet semimajor axis ap‘. 1.2 AU (blue), 1.5 AU (red) and 1.8 AU (orange). Vertical lines denote the value of 
'^p,crit (Eq. 0 for each of the different values of ap. Right: same as left panel, but for a triple configuration as in Fig. In general, LK oscillations are 
entirely suppressed for nip > 2mp crit- For smaller planet mass {nip < ^rit). the eccentricity oscillation amplitude is only slightly modified. 


In our scenario, the abundance of misaligned planets 
around compact binaries depends on the frequency of moderate 
initial inclination stellar triples relative to those with high incli¬ 
nations. High inclination stellar triples may be the progenitors of 
the majority of compact binaries, since the very high eccentrici¬ 
ties reached by the inner binary make orbital decay faster. Our cal¬ 
culations suggest that planets within such high inclination triples 
have less chances of survival during the inner binary’s orbital de¬ 
cay. The efficiency of tidal decay depends on the dissipation time 
scale tv within the stars (see supplementary material). We have 
found that dissipation time scales of order 20 — 50 years can cir¬ 
cularize inner binaries with iin-out,o ^ 78° within a Hubble time, 
but if iin-out,o ~ 70°, then tv — 1 — 5 yr is required. However, 
given the large parameter space in orbital configurations, and the 
uncertainty in realistic values of tv (which may vary during stel¬ 
lar evolution), we cannot discard the possibility that some binaries, 
perhaps still undergoing orbital decay and circularization, may be 
part of moderate-inclination stellar triples, and may therefore be 
candidate hosts to highly misaligned planets. 

An additional caveat to the abundance of misaligned cir- 
cumbinary planets that is not addressed in this work concerns the 
likelihood of planets forming within inclined hierarchical triples 
with flout /flin ~ 100. Planet formation will be limited by disk trun¬ 
cation from inside (at a ~ 3flin) and from outside (at a ~ flout/3) 
|Artymowicz & Lubow| ( |1994] l; [Miranda & Lai] ( |2015^ . Thus, for 
the parameters explored in this paper, planets would be confined 
to form between 1 and 10 AU. In addition to disk truncation and 
warping [Tremaine & Davis] j2014^ , planetesimal dynamics in this 
systems could be affected by the tidal forcing of the inner binary 
and the outer companion, introducing additional complications to 
the formation of planetary cores|Rafikov & Silsbee|(|2015|(;|Silsbee| 
|& RafikovldMSl l. 

As noted before, currently no planets have been detected 
around eclipsing compact (Pin < 5 days) stellar binaries. Our 
work suggests that if planets are able to form within (moderately) 
compact triples, they are likely to survive the tidal shrinkage of 
the central binary, evolving into inclined orbits. The detection of 
these misaligned circumbinary planets may be challenging. The 
planets that survive the orbital decay of the binary lie close to/on 
the Laplace surface, which follows with the precession of the in¬ 
ner binary axis Lin respect to the outer binary axis Lout- The 
coupled precession of the inner binary and the planet orbits will 


produce short-lived “transiting windows”, but these windows ap¬ 
pear periodically over very long time scales [of order 1 /flin-out ~ 
Pin(fflout/oin)^ ~ 10® — 10® years]. An alternative detection strat¬ 
egy is to look for eclipse timing variations. The perturbation on the 
inner binary exerted by a planet of mass nip introduces a timing sig¬ 
nature (on the time scale of the planet’s orbital peri od) of the mag - 
nitude APin ~ (3/87r)Pin(mp/Min)(flin/flp)®^^ Mayer (1990 1 ; 
[Borkovits et al.|p0()3| l. For values of rrip ~ 0.5Mj ~ 0.0005Min, 
Up ~ l.OAU, flin ~ 0.05 AU and Pin ~ 5 days, the maximum 
eclipse timing variation is of order ~ 0.3 s, approaching the noise 
level for some nearby binaries, but in general, still below the detec¬ 
tion limits for most eclipse timing detections [Conroy et al.| ( |^14| ). 
However, short cadence data with current observational capabilities 
might provide enough timing precision to accomplish such mea¬ 
surements. 


We thank Sarah Ballard, Konstantin Batygin, Matthew Hol¬ 
man and Bin Liu for discussions and comments. We also thank the 
referee, Daniel Fabrycky for valuable comments and suggestions. 
This work has been supported in part by NSF grant AST-1211061, 
and NASA grants NNX14AG94G and NNX14AP31G. 

Note added: During of the revision of our manuscript, we 
became aware of a preprint by D. Martin, T. Mazeh and D. Fab¬ 
rycky, which addresses a similar issue (i.e. the dearth of planets 
around compact binaries) as our paper. 
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SUPPLEMENTARY INEOIfMATION 


SI Equations of Motion 

Consider a binary (total mass Min = rno + mi, reduced mass 
Pin = momi/Min and semi-major axis ain) that is a member of 
a hierarchical triple, in which the binary and an outer companion 
of mass Mont orbit each other with a semi-major axis Uout 3^ fflin- 
The shape and orientation of the inner binary orbit are specified 
by the eccentricity vector Oin and the unit vector along the binary’s 
angular momentum direction Lini similarly, the outer companion’s 
orbit is specified by Gout and Lout- A planet orbiting around the 
inner binary has semi-major axis Op, eccentricity vector ep and an¬ 
gular momentum direction Lp. The perturbing potential (per unit 
mass) acting on the planet has contributions from the inner binary, 
{(^in)), and from the outer companion, {{‘f’out)) where the double 
brackets denote time averaging over the orbital periods of the inner 
binary, of the outer companion and of the planet. To quadrupole 
order, these potentials are given by [e.g., [Tremaine et ^(|2009|l] 




^ 1 gMin , _ 2.-3/2 ( /^in ^ f f^in \ 

8 Op V Min J \a-p J 

X 1^1 - 6efn - 3(1 - efn)( Lin • Lp)^ + 15(ein ■ Lp)^j 

(Al) 



Figure Al. Classical Laplace equilibrium for a test particle in a circular 
orbit within a hierarchical stellar triple. The eccentricity of the inner binary 
is set to zero. The inner-to-outer binary inclination fin-out ranges from 0° 
(bottom blue curve) to 80° (top red curve). The dashed portion of the curve 
indicates unstable equilibrium. 


and 


((<E>out)) = - 


1 QMi 
8 a„ 




(- 


\ Q-out 


1^1 - fiCp - 3(1 - ep)( Lout • Lp)^ + 15( Lout • ep)^ 


(A2) 


In our actual calculations, we will set the outer orbit’s eccentric¬ 
ity Bout to zero; this guarantees that higher order octupole terms 
of the potential ((‘bout)) are identically zero. Similarly, setting 
the inner binary to have a mass ratio of unity (which means that 
Uin/Min = 1/4) makes the octupole terms of the inner potential 
((4?in)) vanish. Therefore, the quadrupole potentials given above 
capture the secular dynamics to high accuracy. 

From the potentials ((ibin)) and ((‘bout)), we can derive 
equations of motion for the dimensionless angular momentum vec¬ 
tor of the planet, jp = Lp , and the eccentricity vector Sp 

[see |Tremaine et J)] l |2009[ l; |Liu et aL| ( |2015t ]: 

2 (1 — g2^5/2 ' jp)(jp ^ 

5(6in ' jp)(jp X Gin)j 

“t“ ^ Bout j^(jp ' Lout)(jp X Lout) 5 (Gp • Lout)(Gp X Lout 

(A3a) 




and 

2 (1 — g2)5/2 ^ 

5(0in ' Jp)(®p ^ 

+ 25(ein • jp)" 

— 5(1 — 6in)( Lin • jp) j jp X ep 
Sr." 

“h^^out |^(jp ' Lout)(Gp X Lout) 2jp X Gp 



5(ep • Lout)(jp ^ Lq 


(A3b) 
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Figure A2. Two examples of planet orbital evolution with nearly identical initial conditions showing very different outcomes. As in Fig. 5 in the main text, 
the triple system parameters are ain,o = 0.3 AU, ciout = 30 AU and 2in-out,0 = 83°. The left panels show the eccentricity and inclination evolution of the 
planet (blue) and inner binaiy (green) when ap = 1.6000 AU. The right panels show the same quantities (the evolution of the inner binary is identical in both 
cases) for a planet with ap = 1.6002 AU. The dramatic difference between the two outcomes of the planet’s eccentricity and inclination indicates the chaotic 
nature of planetary orbits undergoing erratic evolution for high-inclination stellar triples. 


where we have introduced the dimensionless quantities 

_ 1 / Min ^ / fflin V 

2 (Min) Up/ ’ 

_ 1 / Mont \ k Op V 

- (1 - (Min) (flont) 

Note that, in the notation of the main body of this article flp-in = 
rtpCin and flp-out — rtpCout- 

By setting djp/dt — dsp/dt = 0 we obtain the equilib¬ 
rium solution for the planet’s orbit under the influence of both in¬ 
ner and outer torques. This is known as “classical Laplace equilib¬ 
rium” in the special case of Cin = 0 |Tremaine et al.| (| 2009| l; [Thmayo| 
|et al.| ( [20T3) l. From Eq. ( |A3a| >, and setting ein = 0, the equation for 
coplanar Laplace equilibrium reads: 

0 = 2ein( Lin ' j)(j ^ Lin) Cout (j ' Lout)(j X Lout) • (A5) 

Using orbital elements, we can write this expression as |Tremaine| 
[eTaLlpOO^ ; [Tamayo et aL| ( [20T^ 

0 — 2ein sin 2 (ip-out iin-out) “t” Cout cos 2ip-out ■ (A6) 


(A4a) 

(A4b) 


We can solve this transcendental equation for ip-out as a func¬ 
tion of the ratio tout/tin, or equivalently, as a function of the 
Laplace radius tl (see Eq. 5 in the main body of this article), since 
{rh/apf = ein/tout. Eig. A1 shows the Laplace surface solution 


as a function of r^/ap for different values of iin-out ranging from 
0° to 90°. This solution of the classical Laplace surface is used in 
Figures 4 and 5 of the main body of the article. 

In our numerical calculations, we directly integrate Eqs. |A3| 
together with the evolution equations of the hierarchical triple as 
the inner binary undergoes LK oscillations with tidal dissipation. 
In principle, this system consists of 24 coupled differential equa¬ 
tions (involving the vectors jin, ein, jp, Op, jout, Sout for the or¬ 
bits, and the spin vectors flo and fli for each of the two central 
stars). We have simplified this system as follows: neglect the evolu¬ 
tion of the outer orbit (valid approximation when the outer angular 
momentum dominates); include only short range forces acting on 
the secondary central star (i.e. the primary star is a non-spinning 
“rigid sphere”); assume psedo-synchronization [e.g., (16)] and or¬ 
bital alignment for the spin vector f2i. This simplification reduces 


the number of equations to 12, while still capturing all the essential 
physics of the problem. 

For the evolution of the vectors ein and Lin = jin/|jin|, 
we use the formalism of|Eggleton & Kiseleva-Eggleton|p()01) and 
[Fabrycky & TremaineH2007y ) 


^jin 

dt 

rfSin 

dt 


rfjin 

dt 

dOin I 

dt I 


in-out 


in-out 


+ 

-f 


djin 

dt 

d©in 

dt 


SRF 

SRF 


(A7a) 

(A7b) 


where the first terms on the right hand side correspond to the con¬ 
servative tidal effect of the tertiary and the second terms to the con¬ 
servative and non-conservative short range forces. In djin /^^IsRF 
and dein /'^^IsRF , an important parameter is the viscous time tv as¬ 
sociated with the dissipation within the stars. This parameter is the 
main source of uncertainty in our calculations. We vary tv between 
~ 1 year to 55 years [the value used by [Fabrycky & Tremainej 
( [2007^ ]. Note that we include dissipation only within one of the 
stars. The effective tidal quality factor Q scales proportionally with 
tv as[Eggleton & Kiseleva-Eggleton|p00T[>: 


Q = 


Qmi tv 


(A8) 


3(l-h2A:)2 i nm ’ 

where we use mi = mo = 0.5Mq, -R*,i = = O.ST?© and k 

(the classical apsidal motion constant) is set to 0.014. We see that 
one to two orders of magnitude of variation in tv is not unreason¬ 
able, given the large degree of uncertainty in the tidal Q values for 
different types of stars. The assumption of pseudo-synchronization 
and spin-orbit alignment also introduces uncertainties, although the 
dominant uncertainty still lies in the value of tv - These uncertain¬ 
ties ultimately affect not the final, circularized orbit of the inner 
binary, but how fast this final state can be reached (see the discus¬ 
sion in the main body of the text). 


S2 Erratic evolution and chaotic behavior 

In Fig. |A2| we show two integrations for planets within the stellr 
triple configuration of Fig. 6 in the main text. The left panels 
show the eccentricity and inclination evolution of a planet with 
Op = 1.6 AU, and the right panels show the same for a planet with 
Up = 1.6002 (about a 0.1% difference in semi-major axis). The 
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Figure A3. Effects of planet mass on the orbital evolution of the inner binary and the planet. Left: Orbital evolution of a binary (oin q = 0.3 AU) and a 
circumbinary planet {a^ = 1.5 AU) in a triple with iin-out ,0 = 73° and Oout = 18 AU where mo + mi = Min = Mout = IMq and rrip = 1 Earth 
mass. The circularization of the inner binary takes place at the same rate as the case with mp = 0. Right: Same as the left panel, but with mp = 1 Saturn 
mass. Circularization is not prevented, but it takes ~ 4 times longer than the rrip = 0 case. 


evolution of both planets is identical until circular orbits become 
unstable, at which point rapid eccentricity growth is triggered, fol¬ 
lowed by highly erratic behavior in the evolution of Cp and ip.out ■ 
The first example “exits” the erratic region and finds a stationary 
state in which Cp and ip-out oscillate around a well-defined value 
with constant amplitude (in this case, the planet lands on a retro¬ 
grade orbit respect to Lout). However, the second example never 
finds “a way out” of the erratic region (which must happen before 
the inner binary circularizes, “locking” the properties of the sys¬ 
tem). Since both these examples go through a phase of extreme ec¬ 
centricities, they are both likely to be ejected by interactions with 
the inner binary, and thus their chances of surviving for a long pe¬ 
riod of time are equally small, regardless of the duration of the 
erratic phase. 


S3 Effect of finite planet mass 


In the calculations above, the mass of the planet rup has been ne¬ 
glected. However, the effect of finite planet mass on the inner bi¬ 
nary can be taken into account in a self-consistent fashion by in¬ 
cluding the tidal potential on the inner binary due to the planet, 
which introduces the additional terms c?j in /df | and dein /df | 

to Eq. |A7| These extra terms have the same functional form as the 
terms arising from the tidal potential of the tertiary djin 
and dein/df|.^^ p Fig. A3 depicts a similar example to that of 


Fig. 3, this time including the effects of nip 7 ^ 0. For a planet of 1 
Farth mass, the results are indistinguishable from the test-particle 
case. However, for a planet of 1 Saturn mass, the effects of the mod¬ 
ified maximum binary eccentricity (see Fig. 7) can be readily seen 
in the efficiency at which the binary shrinks. 


S4 Example of an Al-body integration 

Although long-term direct A-body integrations are computation¬ 
ally costly, one can combine them with the output of the secular 
solutions to study the detailed behavior of the planetary orbit as 
the binary decays. Fig. |A4| shows the A-body result of a dissipa¬ 
tive binary performed using the MERCURY c ode|Chambers|(|1999|>. 
We have added short range forces following [Beauge & Nesvorny| 
( |2012^ , including stellar tides, GR and tidal dissipation. This inte¬ 
gration is started once the binary has undergone appreciable orbital 



10* 10® 10® 
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Figure A4. Example of a dissipative N -body integration of a 4-body sys¬ 
tem. Masses of the binary and companion are mo = mi = 0.5Mq and 
Mout = IMq. The companion semimajor axis is Uout = 18 AU and 
its eccentricity is 0.2. The planet semimajor axis is ap = 1.2 AU and its 
mass is zero. The A-body integration is started once the binary has under¬ 
gone appreciable orbital decay, when ajn = 0.19 AU, Bin = 0.516 and 
*in-out = 73°. The line types are the same as in Fig. 3. 

decay, but before the planet has crossed the Laplace radius vl ; the 
initial condition is obtained using the results of our secular calcula¬ 
tion. In our example, the orbital evolution has been accelerated by 
increasing the tidal dissipation rate. 

As seen from Fig. |A4[ before crossing the Laplace radius, 
the planet’s inclination adiabatically follows that of the inner bi¬ 
nary while maintaining zero eccentricity, in full agreement with our 
secular results. After Up crosses vl, the planet develops a modest 
eccentricity and its inclination decouples from that of the binary. 
Note that since the tertiary companion has a finite eccentricity (0.2) 
in this example, the octupole potential on the planet likely plays an 
important role. This may explain the development of finite planet 
eccentricity for a system with such a modest iin-out- 




























































































































